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The activation cross sections of (d,p), (d,2n), (d, 3n), and (d,2p) reactions on 63,65 Cu were mea- 
sured in the energy range from 4 to 20 MeV using the stacked-foils technique. Then, following 
the available elastic-scattering data analysis that provided the optical potential for reaction cross 
sections calculations, an increased effort has been devoted to the breakup mechanism, the direct re- 
action stripping, and the pre-equilibrium and compound-nucleus cross section calculations, corrected 
for the breakup and stripping decrease of the total reaction cross section. The overall agreement 
between the measured and calculated deuteron activation cross sections proves the correctness of the 
nuclear mechanisms account, next to the simultaneous analysis of the elastic-scattering and reaction 
data. 
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I. INTRODUCTION 

Among the projects of powerful neutron sources for 
nuclear energy generation, the International Fusion Ma- 
terial Irradiation Facility (IFMIF) requests high accu- 
racy deuteron evaluated nuclear data for the assessment 
of induced radioactivity of the accelerator components, 
targets and beam stoppers. The IFMIF facility needs 
such data for estimation of the potential radiation haz- 
ards from the accelerating cavities and beam transport 
elements (Al, Fe, Cr, Cu, Nb) and metal and gaseous 
impurities of the Li loop (Be, C, O, N, Na, K, S, Ca, Fe, 
Cr, Ni) in the energy range from the threshold up to 40 
MeV. However, it is known that the actual experimen- 
tal and evaluated data for deuteron-induced reactions are 
less extensive and accurate than for neutrons, so that fur- 
ther measurements and improved model calculations are 
needed. 

The weak binding energy of the deuteron, £? c /=2.224 



MeV, is responsible for the high complexity of the inter- 
action process that involves a variety of reactions initi- 
ated by the neutrons and protons coming from deuteron 
breakup. Such a wide diversity of nuclear reactions ini- 
tiated by deuteron interaction with nuclei has hampered 
so far the comprehensive analysis involving large A-range 
of targets and incident-energy domain. The difficulties 
to interpret the deuteron-induced reaction data in terms 
of the usual reaction mechanism models have recently 
been investigated [H-@|, looking for a consistent way to 
also include the breakup contribution within the activa- 
tion cross section calculations. Second, the total reaction 
cross sections are less accurately described since, unlike 
the nucleon case, there are no global optical model po- 
tentials (OMP) which describe the scattering data over a 
wide range of nuclei and energies sufficiently well. There- 
fore, a simultaneous analysis of the deuteron elastic scat- 
tering and induced activation [3], which appears es- 
sential for the IFMIF engineering design, is extended in 
the present work for the 63 ' 65 Cu target nuclei. 
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II. MEASUREMENTS 
A. Samples and irradiations 

The variable energy NPI cyclotron provides protons 
and deuterons in energy range 11-37 MeV and 11-20 
MeV, respectively. The irradiation was carried out using 
an external deuteron beam of the NPI cyclotron U-120M 
operating in the negative-ion mode of acceleration. From 
the stripping-foil extractor the beam was delivered to the 
reaction chamber through a beam line consisting of one 
dipole and two quadrupole magnets. 

The incident deuteron energy was determined by com- 
putational procedure based on measured trajectory (the 
frequency and actual extraction radius) of acceleration. 
This procedure was experimentally tested using the acti- 
vation foil method and the surface-barrier-detector tech- 
nique. The energy was determined with a resulting ac- 
curacy of 1.0%, the FWHM spread of the incident beam 
up to 1.8% was observed. 

The activation cross sections induced by deuterons 
bombarding high purity natural aluminum and copper 
foils were measured by a stacked-foil technique. A col- 
limated deuteron beam strikes the stack of foils in a 
Faraday-cup-like reaction chamber enabling the cooling 
of stacked foils without a lost of accuracy in the beam 
current and charge monitoring (10%). To enlarge the 
number of energy bins in the measured excitation func- 
tions and to check the internal consistency of the mea- 
sured data, the foils were stacked with different Al vs Cu 
sequences in two independent runs. The stock of eleven 
Al and eleven Cu foils placed alternately was bombarded 
by a deuteron beam of initial deuteron energy 19.95 MeV, 
with mean beam current 90 nA during an exposure time 
of 15 min. The initial deuteron energy, the beam current 
and the irradiation time for the second run were 20.18 
MeV, 330 nA and 5 min, respectively. The aluminum 
cross section data were reported earlier [U . 

To enable the measurement of the relatively short- 
living isotopes 63 Zn (T 1/2 =38.47 min) and 66 Cu 
(T 1 / 2 =5.12 min), three extra runs were carried out. The 
initial deuteron energy was 19.79, 20.09 and 19.79 MeV, 



TABLE I: Half-lives, main gamma lines and their intensities 
of the measured isotopes 



Isotope T1/2 E 7 I. 







(keV) 


(%) 


yt, Zn 


244.26 d 


1115.5 


50.6 


64 Cu 


12.7 h 


1345.8 


0.473 


62 Zn 


9.186 h 


596.6 


26 






548.4 


15.3 


65 Ni 


2.5172 h 


1481.8 


24 


63 Zn 


38.47 min 


669.6 


8 






962.1 


6.5 


66 Cu 


5.12 min 


1039.2 


9 



the beam current was 0.24, 0.36 and 0.17 /j,A and the 
irradiation time was 6, 5 and 9.3 min, respectively. The 
thickness of the high purity natural Cu and Al foils (pu- 
rity of 99.99%, Goodfcllow product) was 25 and 50 /j,m, 
respectively. Foils were weighted (with a 2% uncertainty) 
to avoid the relatively large uncertainties in the foil thick- 
ness declared by the producer. The mean deuteron en- 
ergy and energy thickness were determined using the 
SRIM 2003 code 0. The overall thickness of the avail- 
able 22 foil stacks covers the excitation-curve range from 
20 to full beam stop. 

In preliminary reports 0, [f| , different initial energies 
were reported due to errors in the orbit calculation of the 
cyclotron operated in the negative-ion mode of accelera- 
tion. In the present report, corrected energy values 19.95 
and 20.18 MeV for the first two runs were established and 
the relevant energies and energy thicknesses of each foil 
were recalculated. 



B. Calculation of cross sections and their errors 

The gamma-rays from the irradiated foils were mea- 
sured repeatedly by two calibrated HPGe detectors of 23 
and 50% efficiency and of FWHM 1.8 keV at 1.3 MeV. 
To provide reliable corrections for the decay, the beam- 
current recorder and 7-ray spectrometer were synchro- 
nized in time. Activated isotopes were identified (see Ta- 
ble HI using nuclear decay data from Ref. [|| . The mea- 
surements with different cooling times lasted up to 100 
days after irradiation. By analyzing the spectra, the re- 
sulting activities at the end of irradiation were obtained. 
The uncertainty of 3% includes statistical errors and the 
uncertainty of the detector-efficiency calibration. 

In the case of short-lived isotope measurements, the 
irradiated Cu foils were immediately measured by the 
HPGe detector with 50% efficiency. To reduce the dead 
time rate caused by the strong annihilation peak ac- 
companying /3 + -decay, the observed Cu foil was situated 
within two iron slides of 1 mm thickness and a lead plate 
of 10 mm thickness was placed between the measured foil 
and the HPGe detector. The detector efficiency was re- 
calibrated according to the experimental conditions using 
a calibrated 152 Eu radioactive source with an uncertainty 
of detector efficiency of 5%. 

The experimental cross sections, given in Table UT1 arc 
shown in Fig. [T] and compared with previously measured 
data [9MT7l| . Their systematic errors are composed of 
current uncertainty (10%), uncertainty of foil thickness 
(2%) and uncertainty of detector efficiency determination 
(2% and 5%, respectively). The mean statistical error 
in activity determination was 2%. Uncertainty of initial 
energy determination was 1%, energy spread of incident 
beam up to 1.8%. Only energy thicknesses are shown in 

Fig.ru 

nat Cu(d, a;) 64 Cu. As the natural copper has two stable 
isotopes, the generation of 64 Cu by irradiation of natu- 
ral copper may proceed via three contribution reactions: 
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TABLE II: Measured reaction cross sections (mb) for deuterons incident on the nat > 63 > 65 Cu nuclei. The mean deuteron energy 
and resolution due to the thickness and straggling of each foil are shown. Statistical uncertainties are given in parentheses for 
the cross sections in units of the last digit. 

Energy Reaction 

(MeV) ""Cu(rf,:r) D ' t Cu DJ Cu(rf, 2n) M Zn ° J Cu(rf, 3n)°^Zn Da Cu(ri, p) DD Cu Da Cu(ri, 2w) DO Zn DO Cu(rf, 2p) 00 Ni 



1.49 (149) 


13.3 (16) 












4.25 


[86) 


85.9 (91) 








9.7 (12) 




4.56 






1.70 (19) 




123 (14) 






5.29 


'751 


195 (23) 








92 (11) 




5.96 


'671 




1.58 (17) 




282 (33) 






6.00 


J66) 








258 (32) 






6.92 


'63) 


214 (24) 








243 (25) 




7.13 


'60) 




8.29 (92) 




322 (36) 






7.69 


'59) 


253 (33) 








445 (53) 




8.21 


'541 




42.0 (48) 




290 (37) 






9.00 


'53) 


216 (28) 








531 (55) 




9.09 


'50) 




103 (12) 




284 (37) 






9.65 


'511 


218 (26) 








675 (81) 




9.91 


'471 




155 (18) 




261 (32) 






10.07 


f47l 




165 (23) 




239 (34) 






10.78 


f47l 


182 (23) 








683 (80) 




11.36 


C451 


188 (25) 








862 (102) 


0.101 (29) 


11.55 


(421 




248 (30) 




207 (24) 






12.25 


(411 




281 (36) 




196 (23) 






12.37 


(431 


161 (17) 








795 (83) 




12.39 


(41) 




269 (52) 




179 (24) 






12.90 


(411 


157 (21) 








880 (104) 


0.200 (42) 


13.69 


(39) 




327 (41) 




165 (20) 






13.82 


(39) 


147 (16) 








940 (95) 


0.355 (73) 


14.31 


(371 


134 (16) 


355 (46) 




166 (25) 


909 (109) 


0.435 (59) 


15.02 


(36) 




353 (83) 




136 (19) 






15.17 


(371 


133 (15) 








949 (95) 


0.502 (64) 


15.62 


(35) 




367 (53) 




138 (19) 






15.63 


(36) 


114 (141 








911 (107) 


0.632 (94) 


16.17 


(341 




380 (47) 










16.44 


(35) 


117 (12) 




0.018 (4) 




901 (91) 


0.750 (83) 


16.87 


(34) 


93 (11) 




0.067 (36) 




851 (100) 


0.70 (11) 


17.37 


(33) 








109 (17) 






17.38 


(32) 




374 (55) 




111 (14) 






17.64 


(33) 


103 (11) 




0.253 (28) 




816 (82) 


0.92(10) 


17.89 


(32) 




382 (54) 




103 (13) 






18.05 


(33) 


91 (11) 




0.645 (89) 




817 (98) 


1.09 (17) 


18.78 


(32) 


102 (12) 




1.44 (15) 




790 (82) 


1.32 (18) 


19.01 


(30) 




366 (66) 




91.6 (122) 






19.17 


(31) 


84 (11) 




2.28 (28) 




730 (88) 


1.27 (18) 


19.49 


(30) 




350 (59) 




89.7 (105) 






19.49 


(29) 




316 (79) 










19.88 


(30) 


97 (12) 




3.72 (38) 




694 (71) 


1.47 (21) 



63 Cu(d,p) (with the threshold E th =0 MeV), 65 Cu(d,t) 
(E t/l =3.767 MeV), and 65 Cu{d,2np) (E th =12.512 MeV). 
Therefore, even for deuteron energies below 20 MeV, we 
have to take into account all contributions and to report 
the measured data to the natural copper as shown in 
Fig. (Ha). 

65 Cu(d,p) 66 Cu is the only possible reaction generating 
66 Cu isotope by deuteron irradiation on natural copper. 
These data are shown in Fig. [lib). 

63 Cu(d, 2n) 63 Zn. There is only one possible way to 
generate 63 Zn in the energy region up to 20 MeV, the 



experimental excitation function is shown in Fig. HJc). 

65 Cu(d, 2n) 65 Zn. The 65 Zn isotope can originate in the 
65 Cu(d, 2n) reaction, and also in the 63 Cu(d, 7) reaction. 
The radioactive capture reaction cross sections are known 
to be very small (of the order of a few hundred /ib or less). 
Hence, the 63 Cu(d, 7) reaction would not be expected to 
contribute appreciably to the measured yield of 65 Zn [sec 
Fig. Hid)]. 

63 Cu(g?, 3n) 62 Zn is the only possible reaction to gener- 
ate 62 Zn in energy region up 20 MeV. The cross sections 
are shown in Fig. [He) ■ 



4 



0.3r 



TABLE III: The parameters of the deuteron optical potential 
for the 63,65 Cu target nuclei. A star used as superscript fol- 
lows the parameters which were changed with respect to the 
optical potential of Daehnick et al. [2J| . 





Potential depths 


Geometry parameters 




(MeV) 


(fin) 


Vj?= 


38.5+0.88Z/A ' -0.26E 


rjj=1.17 






a«=0.734, E<14.5 






=0.709+0.0017E, E>14.5 


w^>= 


-0.014+0.0948E 


r v = 1.325 






a v =0.810 


W* D = 


13.6, E<12 


r D = 1.325 




15.27-0.142E, 12<E<35 


atj=0.770, E<12 




10.3, E>35 


=0.583+0.0156E, 12<E<14.5 






=0.810, E>14.5 


V so = 


=7.33-0.029E 


r so =1.07 






aso=0.66 



65 Cu(e?, 2p) 65 Ni reaction is the only possible way for 
generation of the 65 Ni. Our cross section data (see 
Fig. [TJf)) are the first experimental values except the 
one value around the deuteron energy of 19 MeV. 

Overall, the present experimental data are in satisfac- 
tory agreement with the previous measurements (9Hl7| 
within one standard deviation except for the oldest data 
of Ref . @ and several data of Ref. [l(| • 
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III. ENERGY-DEPENDENT OPTICAL 
POTENTIAL 

The description of the deuteron-nucleus interaction 
represents an important test for both the quality of re- 
action mechanism models and the evaluation of nuclear 
data. The simultaneous analysis of the deuteron elastic 
scattering and induced activation should really involve 
a consistent input of nuclear model calculations 
a prime interest for the optical model potential (OMP) 
parameters being motivated by their further use in the 
analysis of all deuteron interaction cross sections. 

Unfortunately, the few measurements of angular distri- 
butions of elastic scattered deuterons on 63 ' 65 Cu do 
not allow an extended OMP analysis. However, while 
previous OMP analyses on 6 ' 7 LiM [Hj], 27 Al @, HJ, 

54,56,58,natp e ^ 59 Cq and 93 Nb [^Jlhow that no global 

OMP describes sufficiently well the elastic scattering data 
in the energy range up to 20 MeV, but the few pa- 
rameters adjustment (Table IIII[) of the Daehnick et al. 
24, 25] OMP led to a good description of the data for 
the 63 . 65 > nat Cu target nuclei. The comparison of the 
experimental elastic-scattering angular distributions for 
63 ' 65 ' nat Cu [HI and the calculated values obtained by us- 
ing the presently adjusted OMP parameters, the global 
optical potential [2J], and the widely- used TALYS code 
[26| default option based on the Watanabe folding ap- 
proach [27[ (dashed curves) are shown in Fig. At 



FIG. 1: (Color online) Comparison of previous measured data 
[9T447I] and within this work (full circles), the evaluated data 
within the TENDL-2010 library [H] (dotted curves), and cal- 
culated results obtained with TALYS-1.2 code using either 
the whole default input (dash-dotted) or the replacement of 
default deuteron OMP by that of the present work (dashed) , 
for the deuteron interaction with nat . 63 . 65 Cu. 



the same time the measured reaction cross sections for 
deuterons incident on the 63 ' 65 Cu isotopes and natural 
Cu [HI are compared in Fig. [5] with the calculated val- 
ues obtained by using the same potentials as well as the 
evaluated data within the TENDL-2010 library Q]| . One 
may note that the last two calculated excitation func- 
tions underestimate the measured values by at least 20%. 
Finally, the present real-potential diffusibility, imagi- 
nary surface-potential depth and diffusibility, and imagi- 
nary volume-potential depth are compared with the same 
OMP parameters of Daehnick et al. [24j. The elastic- 
scattering cross section calculations have been performed 
using the computer code SCAT2 [29(. 

The particular importance of the deuteron OMP for 
the activation cross section calculations can be seen in 
Fig. [T] through the comparison between the experimental 
data and the calculated results obtained using both the 
present deuteron OMP and the corresponding TALYS 
default option, as well as the TENDL-2010 library [3 
data. There are thus compared the calculated cross sec- 
tions obtained using the default input of TALYS, i.e. the 
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Watanabe folding approach [27j for the deuteron OMP, 
with the results following the replacement of this OMP 
by the parameter set given in Table Mil The differences 
are obvious while it can also be seen that this replace- 
ment still does not lead to a satisfactory description of 
the experimental activation data. Improvements of the 
theoretical analysis by taking into account all reaction 
mechanisms involved in the interaction process are thus 
additionally needed. 
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FIG. 2: (Color online) Comparison of (left) measured 
PJ| deuteron elastic-scattering angular distributions for 
63,65, nat Cu and ca i cmatec i va i ues obtained by using the OMP 
parameters given in Table Mil (solid curves), the global OMP 
[24| (dot-dashed curves), and TALYS 26] default option 
(dashed curves), and (right) measured [28fl reaction cross sec- 
tions for deuterons incident on 63 > 65 .™ a *Cu and calculated val- 
ues obtained by using the same above-mentioned potentials as 
well as the evaluated data within the TENDL-2010 library pH 
(dotted curve). Comparison of particular parameters given in 
Table HIl] (solid curves) and of Daehnick et al. [24] OMP (dot- 
dashed lines) are also shown in the lower part of the right hand 
side. 



IV. DEUTERON BREAKUP 
A. Phenomenological approach 

The interaction of deuterons with the target nuclei pro- 
ceeds largely through direct reaction (DR) processes, for 
deuteron energies below and around the Coulomb bar- 
rier, while with increasing incident energy reaction mech- 
anisms like pre-equilibrium emission (PE) or evapora- 
tion from the fully equilibrated compound nucleus (CN) 
also become important. On the other hand, the breakup 
mechanism is responsible for the enhancement of a large 
variety of reactions along the whole incident-energy range 
and thus its contribution to the activation cross sections 
has to be explicitly taken into account (H-H)]. 

The physical picture of the deuteron breakup in the 
Coulomb and nuclear fields of the target nucleus consid- 
ers two distinct chains, namely the elastic breakup (EB) 
in which the target nucleus remains in its ground state 
and none of the deuteron constituents interacts with it, 
and the inelastic breakup or breakup fusion (BF), where 
one of these deuteron constituents interacts with the tar- 
get nucleus while the remaining one is detected. An 
empirical parametrization of the total proton-emission 
breakup fraction fBu =a Bu/ a R anc ^ the elastic-breakup 
fraction fEB=&EB /<Jr have previously been obtained [2| 
on the basis of experimental systematics (3ll - [35i j . Thus, 
proton-emission spectra and angular distributions from 
deuteron-induced reactions on nuclei from Al to Pb at 
incident energies from 15 to 80 MeV have been studied 
in this respect. However, an energy range of only 15-30 
MeV has been available for the empirical elastic-breakup 
fraction <7eb/vr systematics [33|,[35j]. Their dependence 
on the charge (Z), atomic number (A) of the target nu- 
cleus, and deuteron incident energy (E) was found to be 
0: 

fW n) = 0.087 - 0.0066Z + 0.00163ZA 1/3 + 

0.0017 A 1/3 E - 0.000002Z.E 2 , (1) 



Jeb = 0.031 - 0.0028Z + 0.00051ZA 1/3 + 
0.0005A 1/3 £; - 0.000001ZS 2 . 



(2) 



A comparison with the total proton- and neutron- 
emission breakup cross section parametrization of 
Kalbach [H: 



_f>/« _ K 
"bu ~ 



1 



(13- E) ' ^ d 'P = Kd.n — 18. 

exp- — g — - 



_ (3) 
shows that the former parametrization 2] considers equal 
breakup fractions for proton and neutron emission, but 
also gives all breakup components, i.e. the proton- 
emission breakup total, elastic, and inelastic fraction 

/bj ' —fiiv ' Ieb- The energy dependence of the to- 
tal, elastic, and inelastic proton-emission breakup cross 
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FIG. 3: (Color online) Energy dependence of the total 
(solid curve), elastic (dash-dotted) and inelastic (dot-dashed) 
proton-emission breakup cross sections given by Re f. Bl and 
total proton-emission breakup cross sections of Ref. [3(| (dot- 
ted), for the deuteron interactions with 63 Cu. The corre- 
sponding total reaction cross section is shown by dash-dot- 
dotted curve. 

sections following Ref. Q as well as the total proton- 
emission breakup cross sections |36| for the deuteron in- 
teractions with the 63 ' 65 Cu nuclei are shown in Fig. [3] 
It turns out that, for deuteron incident energies above 
^8 MeV, the predictions for the total proton-emission 
breakup cross sections given by both parametcrizations 
are rather close. However, at the lowest energies the 
total proton-emission breakup cross section provided by 
the latter parametrization [36| become larger than the 
deuteron total reaction cross section. 

Concerning the energy dependence of the inelastic- and 
elastic-breakup components, the interest of the deuteron 
activation cross sections for incident energies up to 60 
MeV has motivated an additional check of the elastic- 
breakup parameterization Q extension beyond the ener- 
gies formerly considered for the derivation of its actual 
form. Actually, as it is shown in Fig. @] for the 63 Cu 
target nucleus, the elastic-breakup cross sections given 
by the empirical parameterization Q decrease with the 
incident energy beyond the energy range within which it 
was established. On the other hand, this trend is opposite 
to that of the total-breakup cross section. Thus, in the 
absence of any available experimental deuteron elastic- 
breakup cross section at incident energies above 30 MeV, 
the correctness of an eventual extrapolation should be 
checked by comparison of the related predictions with 
results of an advanced theory such as the Continuum- 
Discretized Coupled- Channels (CDCC) method [Hil]. 



B. Phenomenological EB versus CDCC formalism 

A detailed description of the CDCC formalism is avail- 
able elsewhere [37H42| | , and hence only a brief description 
of the method is given in the following. 

The breakup component is treated within the CDCC 
formalism as an inelastic excitation of the projectile due 
to the nuclear and Coulomb interactions with the target 



nucleus. Consideration of this excitation is performed 
through the coupling of the projectile unbound excited 
states in the solution of the scattering problem by means 
of the coupled channels approach. The deuteron scatter- 
ing process is analyzed within a three-body model, com- 
prising the two-body projectile and the target, by the 
model Hamiltonian [38j : 

H = K r + K R + V np (r) + U n (R-r/2) + U p (R+r/2). (4) 

Here V v , n is the interaction between the neutron and pro- 
ton [33|, assumed to have a Gaussian shape 

V(r) = -V e~ (r/ro)2 , (5) 

where Vb=72.15 MeV and ro=1.484 fm , were determined 
from the fit of the deuteron binding energy. The vector 
r is the proton-neutron relative coordinate, R is the co- 
ordinate of the center of mass of the p-n pair relative to 
the target nucleus. U p and U n are the proton-target and 
neutron-target interactions, respectively, usually taken 
as the central nuclear part of the proton and neutron 
phenomenological OMPs at half the deuteron incident 
energy, Ed/2- Adjusted Koning-Delaroche [43[ neutron 
and proton global OMPs, in order to obtain a suitable 
description of the deuteron elastic-scattering, have been 
used [44[. The operators K r and Kr are the kinetic en- 
ergies associated with r and R. 

A finite set of coupled equations is obtained by the in- 
troduction of the discretization procedure in which the 
continuum spectrum, truncated at a maximum excita- 
tion energy (S^ax) an d divided into a set of energy (or 
relative momentum) intervals, is represented by a finite 
and discrete set of square-integrable functions. Each bin, 
is represented by a single square-integrable function, cal- 
culated by averaging the scattering states for the p-n rel- 
ative motion within the bin width. Moreover, the p-n 
relative angular momentum is also restricted by consid- 
ering only a limited number of partial waves, in order to 
deal with a finite set of coupled equations. Finally, the 
three-body scattering wave function is expanded over the 
internal states of the deuteron as follows: 

N 

|*(i5)> = (6) 

i=0 

where \<po) is the ground-state wave function and 4>i 
[i 7^ 0) are the averaged (within each bin) continuum 
wave functions. The radial functions Xi(R-) describe the 
projectile-target relative motion for the elastic (i = 0) 
and breakup (i ^ 0) components. Continuum states with 
orbital angular momentum I = 0, 1 and 2 for the p-n rel- 
ative motion were considered. The proton and neutron 
intrinsic spins were ignored for simplicity within calcu- 
lations that were performed with the coupled-channels 
code FRESCO g|. 

The energy dependence of the elastic-breakup cross 
sections provided by the excitation of the continuum 
spectrum, in the case of the deuteron interaction with 
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FIG. 4: (Color online) Energy dependence of the phenomeno- 
logical 2j (dashed curve) and CDCC (solid line) elastic- 
breakup cross sections for deuteron interactions with the 
63 Cu nucleu s I44I1 . The solid circle corresponds to Kleinfeller 
systematics [33j] while there is also shown the total proton- 
emission breakup cross section (dotted curve). 



63 Cu target nucleus, is compared with the prediction of 
empirical systematics Q in Fig. |4] The elastic-breakup 
cross sections corresponding to the Kleinfeller et el. sys- 
tematics (Table 3 of Ref. [33|]) are also shown. The agree- 
ment of the CDCC elastic-breakup cross sections and the 
latter systematics can be considered as validation of the 
present advanced model approach. Moreover, the com- 
parison shown in Fig. |4]points out that the CDCC calcu- 
lations lead to elastic-breakup cross sections which follow 
the total-breakup cross section behavior as well as that 
of the reaction cross section shown in Fig. [3] There- 
fore the present analysis makes clear that the empirical 
parameterization extension for the elastic-breakup cross 
sections beyond the energies considered in this respect 
should be done with caution ^4j. The CDCC method 
thus provides useful initial guidance for the assessment 
of these extrapolation accuracies and may help to im- 
prove the existing phenomenological approaches. How- 
ever, additional experimental deuteron interaction data, 
like elastic-scattering angular distributions and inclusive 
neutron and proton spectra, are needed in order to vali- 
date the parameters involved in the CDCC and complete 
the systematics of the elastic- and total-breakup cross 
sections over enlarged energy and target mass domains. 



C. Inelastic-breakup enhancement of the deuteron 
activation cross sections 

Overall, the deuteron total-breakup cross section 
should be subtracted from the total reaction cross section 
that is shared among different statistical-emission chan- 
nels. On the other hand, the inelastic-breakup processes, 
where one of deuteron constituents interacts with the tar- 
get leading to a secondary composite nucleus, bring con- 
tributions to different reaction channels. The secondary- 
chance emission of particles from the original d-target 
interaction is therefore especially enhanced. Thus, the 



absorbed proton or neutron following the breakup emis- 
sion of a neutron or proton, respectively, contributes to 
the enhancement of the corresponding (d, xn) and (d, xp) 
reaction cross sections. In order to calculate this breakup 
enhancement for, e.g., the (d,xn) reaction cross sections, 
firstly the inelastic-breakup cross sections cr^p were ob- 
tained in the present work by subtracting also the CDCC 
elastic-breakup cross sections from the phenomenologi- 
cal total- breakup cross sections given by Eq. (1). Next, 
they have been multiplied by the ratios crr PtX )/crR 

corresponding to the above-mentioned reactions of the 
absorbed proton with the target nucleus, where ur is the 
proton reaction cross section and x stands for n or 2n 
outgoing channels Q. These ratios have been expressed 
as a function of the deuteron incident energy using the 
Kalbach [46J formula for the center-of-mass system cen- 
troid energy of the deuteron-breakup peak energies of the 
emitted constituents: 



B d T 



l.UZ 



1.5AV3 + 3.1' 



(7) 



In a similar way have been obtained also the inelastic- 
breakup contributions to the (d, p) and (d, 2p) activation 
cross section due to the neutrons absorbed in further in- 
teractions with the target nucleus, i.e. by the (71,7), and 
(n, p) reactions, respectively. The only difference has con- 
sisted in replacing the above-mentioned ratios 0(p. x )l<JR 
by the ratios <J( n , x )/a noni where the non-elastic cross sec- 
tion a non plays the same role for neutrons as <jr for pro- 
tons. 

However, the assumed Gaussian line shape of the 
deuteron-breakup peak energies of the emitted con- 
stituents, that are also showed in Fig. [Sfa) for neutrons, 
have quite large widths. Since the broad approximation 
of this method adopted previously Q for estimation of 
the breakup enhancement, a better estimation is consid- 
ered in the present work. It consists in a convolution of ei- 
ther the ratio <J( n ,x)l a non or crfe^/cR with the Gaussian 
line shape of the deuteron-breakup peak energies of the 
corresponding emitted constituent, for a given deuteron 
incident energy. The cases of deuterons with energies of 
8, 20 and 30 MeV are shown in Fig.[5]Jb) together with the 
cross section ratio <7( Pin )/a\R. There are also shown the 
convolution results at each of these energies, while their 
area corresponds to the inelastic-breakup enhancement 
of the (d,2n) reaction cross sections at these energies. 
These results are more physical, with a realistic incident- 
energy dependence except only for the case involving the 
higher-emission energy side of the Gaussian line shape of 
the deuteron-breakup peak energies. This happens for, 
e.g., deuterons with energies lower than 8 MeV, shown 
for the (d,2n) reaction case in Fig. [5jb) . Since this side 
of the Gaussian line shape could be narrower [46] , differ- 
ent widths for the two halves of the Gaussian distribution 
should be eventually adopted, while otherwise some over- 
estimation may result from using a single width. How- 
ever, the corresponding (d, 2n) reaction cross sections at 
these energies are just above the reaction threshold so 
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FIG. 5: (Color online) (a) The centroid E x of assumed 
Gaussian line shape [4y] for deuteron breakup-peak energies 
of emitted neutrons (solid) and protons (dash-dotted), and 
the corresponding E n ± F/2 values (dashed) calculated for 
deuterons interacting with 63 Cu. (b) The convolution of the 
cross section ratio <T( Pin ) /<?R for the target nucleus 63 Cu (dash- 
dotted curve) with the Gaussian line shape of the deuteron- 
breakup peak energies of the corresponding emitted protons, 
for deuterons with energies of 8, 20 and 30 MeV (solid curves, 
with the incident energy noted above their maxima), and the 
convolution results at each deuteron energy (dashed curves). 



that we have taken into account the inelastic-breakup 
enhancement of the {d, 2n) reaction cross sections only 
above these energies. 



V. ONE-NUCLEON TRANSFER REACTIONS 

Apart from the breakup contributions to deuteron in- 
teractions, the direct reaction mechanisms like stripping 
and pick-up have to be properly considered in order to 
describe the low energy side of the (d, p), (d, n) and (d, t) 
excitation functions 23]. In the present work the DR 
contribution to the (d, p) and (d, t) reaction reaction cross 
sections, through population of the low-lying discrete 
levels of 64,66 Cu residual nuclei, was calculated using 
the code FRESCO [H[ based on the Coupled-Reaction 
Channels (CRC) method. The post form distorted-wave 
transition amplitude with finite-range interaction has 
been chosen. 



A. The one-nucleon stripping 

The (rf, p) reaction has been in a large extent of critical 
importance for the study of nuclear structure. Actually, 
the spectroscopic factors extracted from the analysis of 
experimental angular distributions of the corresponding 
emitted protons did contribute to the validation of the 
nuclear shell model, considering that the neutron from 
the deuteron is transferred to a single-particle state of 
the residual nucleus. Consequently, the rich systematics 
of the achieved experimental spectroscopic factors makes 
possible the calculation of almost total stripping cross- 
section contribution to the deuteron activation cross sec- 
tion. 

The above-mentioned deuteron phenomcnological 
OMP parameter set (Table IIII[) has been used for the 
incident channel, while the Koning-Delaroche [43| OMP 
global parameters have been used for protons interac- 
tions with the residual nuclei. The n-p interaction in the 
deuteron has been described with the potential V np given 
by Eq. (O and the neutron bound states were generated 
in a Woods-Saxon real potential with the global values 
of a reduced radius of 1.25 fm and diffuseness of 0.65 
fm, while its depth has been adjusted to reproduce the 
nucleon binding energies in the residual nuclei. 

The present calculations of the single-neutron stripping 
(d,p) reaction cross section have involved transitions to 
104 final states of the odd-odd residual nucleus 64 Cu and 
to 81 final states of the similar residual nucleus 66 Cu. The 
spectroscopic factors that were obtained experimentally 
from proton angular-distribution measurements for these 
states up to ~5 MeV, as given in Table II of Ref. [47j j . 
Table I of Ref. [48| , and [43 . |50| , have been considered in 
this respect. One may note that a lower number of final 
states, i.e. 63 for 64 Cu and 52 for 66 Cu, extending up 
to ^3 MeV, were used within a preliminary stage of this 
work Q , which made necessary an additional assumption 
concerning the DR contribution from the states at higher 
energies. The present increase of the final states taken 
into account for the DR contribution to the (d, p) reaction 
cross sections makes a similar assumption no longer nec- 
essary. In spite of the corresponding results being rather 
close to the former ones Q , an increased accuracy is now 
obtained for this activation component, even better than 
-5%. 

The resulted stripping components of the {d,p) reac- 
tion excitation functions are essential for the description 
of the experimental data shown in Fig. [BJ This statement 
is valid also concerning the maxima of these excitation 
functions at deuteron incident energies £?d~8 MeV. 

B. The one-nucleon pick-up 

The excitation function corresponding to 64 Cu nucleus 
production from deuterons incident on a natural cop- 
per target includes contributions from both stripping 
63 Cu(d,p) and pick-up 65 Cu(d, t) direct reaction mech- 
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anisms. Actually, the lowest energy side of a (d, t) exci- 
tation function, between its threshold and those of the 
(d,dn) and (d,p2n) reactions, can be described exclu- 
sively by the pick-up reactions as it is shown, e.g., in Fig. 
3 of Ref. HI. 

Therefore, we have used in the present work the above- 
mentioned deuteron phenomenological OMP to describe 
the incident channel of the (d, t) reaction, in a similar 
way to the stripping calculation, while the Becchetti and 
Greenlees pH. l5l| OMP has been used for the emitted 
tritons. 

The d-n interaction in triton has been described with 
a 3 He potential 52] of Woods-Saxon shape: 



V(r) = -V 



1 



1 + e ( r ~ r °y a ' 



(8) 



where the parameters Vo=77.71 MeV, ro=1.008 fm and 
a=0.75 fm were determined by a fit of the 5.50 MeV 3 He 
binding energy (relatively close to 6.3 MeV corresponding 
to 3 H). The transfered-neutron bound states were also 
generated in a Woods-Saxon real potential, with global 
reduced radius of 1.25 fm, diffuseness of 0.65 fm, and 
the depth adjusted to the nucleon binding energies in 
the residual nuclei. The experimental spectroscopic fac- 
tors obtained by analysis of triton angular distributions 
related to the population of 14 discrete levels of the resid- 
ual nucleus 64 Cu [49[ have been used for calculation of 
the pick-up transition amplitude. 

The contribution thus obtained of the 65 Cu(d, t) re- 
action to the Cu production excitation function, cor- 
rected for the isotopic ratio 65 Cu/ 63 Cu, is shown in 
Fig- EUa). As expected [23| . the (d, t) activation cross 
sections is the essential contribution among the processes 
induced by deuterons on the 65 Cu isotope, within the 
natural copper target, at incident energies between 7-20 
MeV. The statistical emission through the (d, 2np) and 
(d,nd) reactions, also shown in Fig. |6](a), become signif- 
icant at higher energies. 



VI. STATISTICAL PARTICLE EMISSION 

The reaction mechanisms such as pre-equilibrium emis- 
sion (PE) or evaporation from the fully equilibrated com- 
pound nucleus (CN), become important when the in- 
cident energy is increased above the Coulomb barrier, 
when the interaction of deuterons with the target nu- 
clei proceeds largely by DR processes. The related cross 
sections have been analyzed in this work by using the de- 
fault model parameters (except for the deuteron OMP in 
Table IIII[) of the widely-used computer code TALYS as 
well as a local consistent parameter set developed in cal- 
culations with the PE+CN code STAPRE-H |H3| taking 
into account also the breakup and DR results discussed 
above. 

The local analysis results obviously have a higher accu- 
racy while the global predictions may be useful for an un- 
derstanding of unexpected differences between measured 



and calculated cross sections. The main assumptions and 
parameters involved in this work for the sets of global and 
local calculations have recently been described elsewhere 
[55|, only some points specific to the mass range A>60 are 
given here. A further note should concern the fact that 
similar input parameter sets and calculations have been 
used to obtain the breakup-enhancement due to one of 
deuteron constituents interacting with the target nucleus 
and leading to a secondary composite nucleus, with final 
contributions to different reaction channels as discussed 
in Sec. [TYl 

The deuteron phenomenological optical model param- 
eter set given in Table Hill has been used for the incident 
channel. The nucleon optical potential of Koning and 
Delaroche [HI, used by default in the TALYS code, has 
obviously been the first option. However, a basic point 
revealed by these authors is that their global potential 
does not reproduce the minimum around the neutron en- 
ergy of 1-2 MeV for the total neutron cross sections of 
the A^60 nuclei. Following also their comment on the 
constant geometry parameters which may be responsible 
for this aspect, we have applied the SPRT method [HiJ 
for determination of the OMP parameters over a wide 
neutron energy range through analysis of the s- and p- 
wave neutron strength functions, the potential scattering 
radius R' and the energy dependence of the total cross 
section ax{E). The recent RIPL-3 recommendations (25| 
for the low-energy neutron scattering properties and the 
available measured <tt data have been used in this re- 
spect, and we found that it is necessary to consider the 
energy dependence of the real potential geometry at lower 
energies given in Ref. 



55|. 



These potentials were also used for the calculation of 
the collective inelastic scattering cross sections by means 
of the direct-interaction distorted-wave Born approxi- 
mation (DWBA) method and a local version of the com- 
puter code DWUCK4 [56[ . The weak coupling model was 
adopted for the odd nuclei 55 Mn and 63,65 Cu using the 
collective state parameters of Kalbach 57J. Typical ra- 
tios of the direct inelastic scattering to the total reaction 
cross sections in the energy range from few to 60 MeV 
decrease from ~11 to 5%, for the 55 Mn nucleus, and from 
^8 to 3% for the Cu isotopes. 

The OMP of Koning and Delaroche was also consid- 
ered for the calculation of proton transmission coefficients 
on the residual nuclei, i.e. the isotopes of Cu and Ni, 
while a previous trial of this potential concerned the pro- 
ton reaction cross sections an [H^ ]. Actually our local 
analysis involved the isotopes of Mn, Fe, Co, Ni, Cu and 
Zn, for lower energies important in statistical emission 
from excited nuclei. In order to obtain the agreement 
with the corresponding <7r data we have found it neces- 
sary to replace the constant real-potential diffusivity by 
the energy-dependent form ay—0- 463+0. 01E up to 20 
MeV for 58 Ni, where the energy E is in MeV and the dif- 
fusivity is in fm. A final validation of both the original 
OMP and the additional energy-dependent ay has been 
obtained by analysis of the available (p, 7) and (p, n) re- 
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action data up to E p ~12 MeV on Ni isotopes while the 
other statistical model parameters are the same as in the 
rest of the present work. 

The optical potential which is used in this work for cal- 
culation of the a-particle transmission coefficients was es- 
tablished previously 59] for emitted a-particles, and sup- 
ported recently by semi-microscopic analysis for ^4~90 
nuclei |60j. On the other hand, by comparison of the 
present calculations and measured data [30J for the tar- 
get nuclei 63 ' 65 Cu we found that the real well diffusc- 
ness an of the above-mentioned global OMP should be 
changed to 0.67 fm. This reduction is rather similar 
to that found necessary for the target nuclei 59 Co, and 
58 ' 60 ' 62 Ni [6l|. Moreover, a final remark concerns the fact 
that the same OMP parameter sets were employed within 
both the PE generalized [53[ Geometry-Dependent Hy- 
brid (GDH) model (62}, and the CN statistical model. 

The nuclear level densities were derived on the basis 
of the back-shifted Fermi gas (BSFG) formula (63|, for 
the excitation energies below the neutron-binding energy, 
with small adjustments of the parameters a and A [64] 
obtained by a fit of more recent experimental low-lying 
discrete levels [65J and s-wave nucleon resonance spacings 
Do [13 ■ Above the neutron binding we took into account 
the washing out of shell effects within the approach of 
Ignatyuk et al. (6(| and Junghans et al. |s3|j an d using 
the method of Koning and Chadwick [68[ for fixing the 
appropriate shell correction energy. A transition range 
from the BSFG formula description to the higher energy 
approach has been chosen between the neutron binding 
energy and the excitation energy of 15 MeV, mainly in or- 
der to have a smooth connection. On the other hand, the 
spin distribution has been determined by a variable ratio 
I /I r of the nuclear moment of inertia to its rigid-body 
value, between 0.5 for ground states, 0.75 at the neu- 
tron binding energy, and 1 around the excitation energy 
of 15 MeV. Concerning the particle- hole state density 
which for the PE description plays the same role as the 
nuclear-level densit y fo r statistical model calculations, a 
composite formula [69| was used within the GDH model 
with no free parameters except for the a-particle state 
density ,g Q =A/10.36 MeV" 1 [zO]. 

The modified energy-dependent Breit-Wigner 
(EDBW) model [zE [z3 was used for the electric dipole 
7-ray strength functions /ei(2?~) of main importance 
for calculation of the 7-ray transmission coefficients. 
The corresponding /ei(.E7 7 ) values have been checked 
within the calculations of capture cross sections of Mn 
and Cu isotopes in the neutron energy range from keV 
to 3-4 MeV, by using the OMP and nuclear level densit 
parameters described above and global estimations [7T 
of the 7-ray strength functions for multipoles A<3. 

Thus we found that the /ei(-E 7 ) strength functions 

corresponding to the experimental [25[ average radiative 
widths r^g P provide an accurate description of the 
capture data for the Cu isotopes. Finally, the accuracy 
of the 7-ray strength functions adopted in this work is 
also shown by the above-mentioned analysis of the (p, 7) 



reaction cross sections. 

Formally, no free parameter is involved for the 
PE description within the corresponding generalized 
GDH model except for a-particle emission, the above- 
mentioned s. p. 1. -density and the pre-formation proba- 
bility ip [z3] with a value of 0.2 used in the present 
work. However, a particular comment concerns the ini- 
tial configuration of excited particles (p) and holes (h), 
for deuteron-induced reactions in the present case. Simi- 
lar careful studies [U [H, [H, [74[ pointed out that 3p-lh 
or 2p-lh may be a suitable choice for this configuration. 
Our calculations show that the latter one gives the best 
agreement between the measured and calculated reaction 
cross sections. 



VII. RESULTS AND DISCUSSION 

The comparison of the measured and calculated (d, p) 
reaction cross sections of raa *. 63 > 65 Cu is shown in Fig. [51 
including the present experimental data already shown 
in Fig. [TJ and the global and local analysis results. For 
the local analysis both components of the final activation 
are shown, i.e. the DR cross sections provided by the 
code FRESCO and the PE+CN contributions supplied 
by STAPRE-H. The latter is alone rather close to the 
TALYS predictions. The local approach has led to much 
better agreement with the present (d, p) reaction data 
0, [|[ especially due to the stripping DR contribution. 

In order to obtain a complete description of the (d, 2p) 
reaction cross sections, we have started by taking into 
account also the neutrons which, following the breakup 
proton emission, are absorbed in further interactions with 
the target nucleus. The cross section a p BF 22[ has been 
considered in this respect (Fig. [5]) as well as the corre- 
sponding fraction leading to the above-mentioned reac- 
tions. These fractions have been obtained by using the 
ratios of the most recently evaluated [zH {n,p) and (n, a) 
reaction cross sections, respectively, to the neutron reac- 
tion cross sections provided by the neutron global OMP 
[43l | . A similar procedure has been followed in order to 
obtain the contribution to the (d, 2n) and (d, 3n) reaction 
cross sections due to the protons which, following the 
breakup neutron emission, are absorbed in further inter- 
actions with the target nucleus and described by the cross 
section <J%p. The only difference in this case concerns the 
(p, n) reaction cross sections in the incident energy range 
up to 30 MeV, which have been obtained by PE+CN cal- 
culations using the computer code STAPRE-H and the 
consistent local parameter set described above. All in- 
termediary and ultimate reaction cross sections shown 
in Fig. [5][c-e) indicate that they may contribute up to 
50% of the activation cross sections for deuteron incident 
energies of ~25 MeV. 

The contribution due to the breakup proton, added 
to the PE+CN components provided by STAPRE-H, 
describe rather well the measured cross sections of the 
(d,2p) reaction as shown in Fig. [BJf). Similarly, the 
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FIG. 6: (Color online) Comparison of measured data already shown in Fig. Q] and present analysis results (solid) taking into 
account the deuteron inelastic breakup (dashed), the DR (dash-dotted for (d,p) reactions and short dash-dotted for (d,t) 
reaction), and the PE+CN (dotted and short-dotted for (d, 2np) reaction) mechanism contributions to the deuteron interaction 
with ™ at > 63 . 65 Cu target nuclei. 



breakup neutron emission plays the same role for the 
(d, 2n) reaction as shown in Fig. [UJc,d). Their weight is 
obviously increasing with the incident energy since all re- 
actions involved, following the deuteron breakup, within 
the second step of these processes have negative Q- values. 

Finally, all activation data of deuteron-induced reac- 
tions on 63,65 Cu have been properly described, making 
obvious the usefulness of the concurrent description of 



all reaction channels as well as the simultaneous analysis 
of the deuteron elastic scattering and induced activation. 
A particular underprediction has concerned however the 
65 Cu(<i,p) 66 Cu reaction cross sections, only their energy 
dependence being well described. A first comment may 
concern in this respect the fact that, although the TALYS 
and TENDL calculations do include a breakup compo- 
nent in all (d,n) and {d, p) reaction channels, the sys- 
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tematical relations for its strength does not show enough 
predictive power in this particular case. On the other 
hand, also the lower number of the known spectroscopic 
factors corresponding to the discrete states of the odd- 
odd residual nucleus 66 Cu, taken into account for the DR 
contributions, may explain this underprediction. 

VIII. SUMMARY 

The cross section values for deuteron-induced reac- 
tions on natural Cu were determined for the reac- 
tions ™ at Cu(d,a;) 64 Cu, 65 Cu(d,p) 66 Cu, 63 Cu(d, 2n) 63 Zn, 
65 Cu(d, 2n) 65 Zn, 63 Cu(d, 3n) 62 Zn and 65 Cu(d, 2p) 65 Ni at 
deuteron energies up to 20 MeV. Resulting cross section 
data are in good agreement with the major part of pre- 
vious reported experiments. 

Following a previous extended analysis of elastic- 
scattering, breakup and direct-reactions of deuterons on 
63 ' 65 Cu, for energies from 3 to 60 MeV 0, the pre- 
equilibrium and statistical emissions have been consid- 
ered in the same energy range. The related cross sec- 
tions have been analyzed by using the default model pa- 
rameters (except for the deuteron OMP in Table UTTj) of 
the widely-used computer code TALYS as well as a local 
consistent parameter set developed in calculations with 
the PE+CN code STAPRE-H taking into account also 
the breakup and DR results formerly discussed. The lo- 
cal approach has led to much better agreement with the 
present (d,p) reaction data especially due to the model 



calculation of the stripping DR contribution. 

Consideration of the deuteron breakup plays a key role 
for the reaction channels adding a second emitted parti- 
cle to the first one. Thus, in order to obtain a complete 
description of the (d,p), (<i, 2n), (d, 3n), and (d, 2p) reac- 
tion cross sections, we have taken into account also the 
neutrons which, following the breakup proton emission, 
are absorbed in further interactions with the target nu- 
cleus. Finally, all deuteron-induced reactions on 63 ' 65 Cu, 
including the present data measured at 20 MeV deuteron 
energy, have been properly described due to a simultane- 
ous analysis of the elastic-scattering and reaction data. 
A similar analysis will be further considered for a sys- 
tematical evaluation of the deuteron activation of other 
medium-mass nuclei. 
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